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Abstract. LS I +61 303 is one of the most studied X-ray binary systems because of its two 
peculiarities: On the one hand being the probable counterpart of the variable gamma ray source 
2CG 135+01 0112] an d on the other hand being a periodic radio source |3]. The recent discovery 
of a radio emitting jet extending ca. 200 AU at both sides of a central core 01 in all evidence has 
shown the occurrence of accretion/ejection processes in this system. However, the radio outbursts 
do not occur at periastron passage, where the accretion is at its maximum, but several days later. 
In addition, when the gamma-ray emission of 2CG 135+01 is examined along the orbital phase of 
LS I +61 303 one sees that this emission seems to peak at periastron passage [5]. Here in detail we 
analyse the trend of the gamma-ray data versus orbital phase and discuss the delay between peaks 
at gamma-rays and in the radio band within the framework of a two-peak accretion/ejection model 
proposed by Taylor et al. |6] and further developed by Marti & Paredes |7]. 
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X THE BINARY STELLAR SYSTEM 

a' 

The binary stellar system LS I +61 303 has an orbital period of 26.4960 days and a rather 
eccentric orbit (e=0.7) j l, l3il7ll9i Holl . The visible companion is an early type, rapidly 
rotating BOVe star ISlllPj]. Be stars have a dense and slow (<100 km s _1 ) disk-like wind 
with a wind density distribution following a power law; a long monitoring of the H a line 
emission has indicated that in the system LS I +61 303 are present variations in the mass 
loss with a period of about 4.6 years JTl, 12, 3- 



Through this dense, structured and variable envelope the compact object travels and 
accretes (see sketch in Fig. 1). It is still an open issue whether the compact object in 
this system is a neutron star or a black hole. Recent optical observations [10] give 
a mass function / in the range 0.003 < f < 0.027. The maximum value of / for an 
inclination i = 38° and a mass of the Be star M = 18M implies a mass of the compact 
object Mx < 3.8 M ; since the assumed upper limit for a stable neutron star is 3M Q , the 
presence of a black hole in the system cannot be ruled out 10] . 




FIGURE 1. Sketch of the binary system LS I +61 303: Orbital phases of the compact object travelling 
through the wind of the Be companion (see text). 



RADIO EMISSION 

One of the most unusual aspects of the radio emission is that it exhibits two periodicities: 
a 26.4960 day periodic nonthermal outburst and a 4.6 years modulation of the outburst 
peak flux IBllal. The last period is clearly correlated with the mass loss of the Be star 
ill ll H2L Il3h whereas the shorter periodicity corresponds to the orbital period of the 
binary system JEHU- The variations of the Be equatorial disk influence, besides the 
outburst peak flux, the orbital phase where the outbursts occur and give rise to the 
observed broad distribution r adio-outbursts=O. 45-0.95 [ 14] where the phase is calculated 
for t = JD 2 443 366.775 and P = 26.4960 days §]. 

As shown in Fig. 1 the <j) at periastron passage is 0.2 1 1011 - Therefore one of the 
fundamental questions concerning the periodic radio outbursts of LS 1+61 303 has been: 
If the maximum accretion must occur where the density is highest, why are the radio 
outbursts delayed with respect to the periastron passage? 

On the other hand the radio emission has been resolved with VLBI, EVN and MER- 
LIN lfl5L \ldi 0, The morphology of the radio emission in Fig. 2a (at 0=0.7) is 
typical for a microquasar, i.e. jet-like, and constitutes the observational evidence of the 
occurrence of accretion/ejection processes in this system. Moreover, the morphology in 
Fig. 2a is S-shaped and quite similar to that of the precessing jet of SS 433 lfl8ll . The 
receding jet is attenuated by Doppler de-boosting but still above the noise limit of the 
image. One day later (Fig. 2b) the precession suggested in the first MERLIN image 
becomes evident: a new feature is present oriented to the North-East at a position an- 
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FIGURE 2. a) MERLIN image of LS I +61 303 at 5 GHz obtained on 2001 April 22. North is up and 
East is to the left. The synthesized beam has a size of 5 1 x 58 mas, with a PA of 17°. The contour levels are 
at -3, 3, 4, 5, 6, 7, 8, 9, 10, 20, 40, 80, and 160(7, being d=0.14 mJy beam" 1 . The S-shaped morphology 
strongly recalls the precessing jet of SS 433, whose simulated radio emission (Fig. 6b in Hjellming & 
Johnston fljio is given in the small box. b) Same as before but for the April 23 run. The synthesized beam 
has a size of 39 x 49 mas, with a PA of — 10°. The contour levels are the same as those used in the April 
22 image but up to 320a, with ct=0.12 mJy beam -1 fl. 



gle (PA) of 67°; the counter-jet is severely de-boosted. The Northwest-Southeast jet of 
Fig. 2a has a PA=124°. Therefore a quite large rotation has occurred in only 24 hours. 



GAMMA-RAY EMISSION 

In 1978 Gregory and Taylor Jll reported about the discovery of a radio source (later 
on associated to LS 1 +61 303) within the lo error circle of the COS B y-ray source 
2CG 135+01. This association remained however controversial because of the presence 
of the quasar QSO 0241+622 within the relatively large COS B error box. The position 
of this gamma-ray source given as 2EG J0241+6119 in the Second EGRET Catalog 
is / = 135°. 58, b = 1°.13; the radius of the 95% confidence error contour of about 
13 arcminutes has finally ruled out the possible identification with QSO 0241+622, 64 
arcminutes away. The position is, on the contrary compatible with LS I +61 303 only 
8 arcminutes away II 1911 . In 1998 Tavani and collaborators establish the possibility of 
variability of 2CG 135+01 on timescales of days [2]. Massi [5] examines the EGRET 
data as a function of the orbital phase and notices the clustering of high flux around the 
<j) interval 0.2-0.5 estimated to be the periastron passage. 

Very recently Casares and collaborators 111 Oil have finally confirmed periastron pas- 
sage to be at 0=0.2. Let us, therefore examine the plot of Fig. 3. The plot begins at the es- 
tablished periastron passage =0.2 and shows the follow-up of the EGRET gamma-ray 
emission along one full orbit. At epoch JD 2 450 334 (i.e. circles in the plot, with empty 




circles indicating upper limits) the orbit has been well sampled at all phases: A clear 
peak is centered at periastron passage 0.2 and 1.2. At a previous epoch (JD 2449045; 
triangles in the plot) the sampling is incomplete, still the data show an increase near 
periastron passage at ~0.3, and a peak at ~0.5. The 3 squares refer to a third epoch 
(JD 2449471). 



THE TWO-PEAK ACCRETION/EJECTION MODEL 

Taylor et al. |0] and Marti & Paredes 0] have modelled the properties of this system 
in terms of an accretion rate M oc (where p w i n d is the density of the Be star wind 

V rel 

and v re i is the relative speed between the accretor and the wind) which develops two 
peaks for high eccentricities: the highest peak corresponds to the periastron passage 
(highest density), while the second peak occurs when the drop in the relative velocity 
v re i compensates (because of the inverse cube dependence) the decrease in density. Marti 
& Paredes [7] have shown that both peaks are above the Eddington limit and therefore 
one expects that matter is ejected twice within the 26.5 days interval. During the first 
ejection, due to the proximity to the Be star, no radio emission but high energy emission 
is expected because of severe inverse Compton losses. 

As a matter of fact whereas the few gamma-ray data available show high levels 



TABLE 1. Luminosity (erg/s) Q 



L r (>100MeV) L x Ladio 

3.1 xlO 35 (l-6)x!0 34 (1-I7)xl0 31 



of gamma-ray emission around periastron passage (previous section), radio outbursts 
were indeed never observed at periastron passage in more than 20 years of radio flux 
measurements. Moreover, Bosch-Ramon and Paredes [20] have proposed a numerical 
model based on inverse Compton scattering where the relativistic electrons in the jet 
are exposed to stellar photons (external Compton) as well as to synchrotron photons 
(synchrotron self Compton). This model considers accretion variations along the orbit 
and predicts indeed a gamma-ray peak at periastron passage where the accretion is 
higher. 

At the second accretion peak the compact object is enough far away from the Be 
star, so that the inverse Compton losses are small and electrons can propagate out of the 
orbital plane. Then an expanding double radio source should be observed, what in fact 
has been observed (Fig. 1) by MERLIN at orbital phase 0.7 (apoastron). Interesting in 
this respect is the gamma-ray peak at ~ 0.5 that could originate from a second ejection 
which occurred still close enough to the Be star. 

CONCLUSIONS 

It is clear that new gamma-ray observations are desirable to confirm the periodical 
gamma-ray flares at periastron passage predicted by the model. For a second ejection 
still close enough to the Be star in order to generate a second gamma-ray flare, joint ra- 
dio observations would allow to determine important physical parameters. The radio and 
gamma-ray peaks are expected to be delayed: When the electrons start to produce syn- 
chrotron emission in the radio band, this emission is absorbed by the thermal electrons 
of the Be wind. Such a delay strongly depends on the distance from the Be star. 

In ejections distant enough from the Be star it could be verified if also disk photons 
are upscattered by the relativistic electrons and contribute to the gamma-ray emission 
iEHi . such a contribution would finally explain why LS 1+61 303 is subluminous in the 
X-ray range (see Table 1, adapted from J22J]). 

Moreover, on the basis of the precession shown in Fig. 2, we suggest that the fast 
precession bringing the jet intermittently closer and farther to the line of sight should 
produce noticeable variable y-ray emission. In fact, the amplification due to the Doppler 
factor for Compton scattering of stellar photons b y th e relativistic electrons of the jet 
is even higher than that for synchrotron emission l!2lL l4ll. In conclusion LS I +61 303 
becomes the ideal laboratory to test the recently proposed model for microblazars with 
INTEGRAL and MERLIN observations now and by AGILE and GLAST in the future. 
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